Quantum fidelity for degenerate groundstates in quantum phase transitions 
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Spontaneous symmetry breaking mechanism in quantum phase transitions manifests the existence of degen- 
erate groundstates in broken symmetry phases. To detect such degenerate groundstates, we introduce a quantum 
fidelity as an overlap measurement between system groundstates and an arbitrary reference state. This quantum 
fidelity is shown a multiple bifurcation as an indicator of quantum phase transitions, without knowing any de- 
tailed broken symmetry, between a broken symmetry phase and symmetry phases as well as between a broken 
symmetry phase and other broken symmetry phases, when a system parameter crosses its critical value (i.e., 
multiple bifurcation point). Each order parameter, characterizing a broken symmetry phase, from each of de- 
generate groundstates is shown similar multiple bifurcation behavior. Furthermore, to complete the description 
of an ordered phase, it is possible to specify how each order parameter from each of degenerate groundstates 
transforms under a symmetry group that is possessed by the Hamiltonian because each order parameter is in- 
variant under only a subgroup of the symmetry group although the Hamiltonian remains invariant under the full 
symmetry group. Examples are given in the quantum g-state Potts models with a transverse magnetic field by 
employing the tensor network algorithms based on infinite-size lattices. For any q, a general relation between 
the local order parameters is found to clearly show the subgroup of the X q symmetry group. In addition, we 
systematically discuss the criticality in the g-state Potts model. 

PACS numbers: 03.67.-a, 05.30.Rt, 05.50.+q, 75.40.Cx 



I. INTRODUCTION 

Quantum phase transitions (QPTs) have attracted much 
attention to understand a relationship with quantum 
information 12 . Compared to local order parameters in 
the conventional Landau-Ginzburg-Wilson paradigm based 
on the spontaneous symmetry breaking mechanism, in con- 
nections between quantum information and QPTs, quantum 
entanglement, i.e., a purely quantum correlation being absent 
in classical systems, can be used as an indicator of QPTs 
driven by quantum fluctuations in quantum many-body 
systems 3 . Quantum fidelity, based on the basic notions of 
quantum mechanics on quantum measurement, has also 
provided an another way to characterize QPTs 4 7 . In the 
last few years, various quantum fidelity approaches such 
as fidelity per lattice site (FLS) 4 , reduced fidelity 8 , fidelity 
susceptibility 9 , density-functional fidelity 9 , and operator 
fidelity-^, have been suggested and implemented to explore 
QPTs. 

Actually, the fidelity is a measure of similarity between two 
quantum states by defining a overlap function between them. 
The fact that groundstates in different phases should be or- 
thogonal due to their distinguishability in the thermodynamic 
limit allows a fidelity between quantum many-body states in 
different phases signaling QPTs because an abrupt change of 
the fidelity is expected across a critical point in the thermo- 
dynamic limit 4 - - — . Thus, the fidelity has great advantages 
to characterize the QPTs in a variety of quantum lattice sys- 
tems because the groundstate of a system undergoes a drastic 
change in its structure at a critical point, regardless of of what 
type of internal orders are present in quantum many-body 
states. Especially, the groundstate FLS has been manifested 
to capture drastic changes of the groundstate wave functions 
around a critical point even for those cannot be described in 
the framework of Landau-Ginzburg-Wilson theory, such as a 
Beresinskii-Kosterlitz-Thouless (BKT) transtion and topolog- 



ical QPTs in quantum lattice many -body systems 11 . 

Even though such latest advances in understanding QPTs 
have been made significantly, however, understanding directly 
degenerate groundstates originated from a spontaneous sym- 
metry breaking as a heart of the Landau-Ginzburg-Wilson the- 
ory has still been remained in unexplored research regimes. 
Also, recently developed tensor network (TN) algorithms^— 
in quantum lattice systems have made it possible to explore di- 
rectly degenerate groundstates with a randomly chosen initial 
state subject to an imaginary time evolution. Based on the ten- 
sor network algorithm, indeed, Zhao et al^, for the first time, 
have detected a doubly degenerate groundstate by means of 
FLS bifurcations in quantum Ising model and spin- 1/2 XYX 
model with transverse magnetic field. Also, in various spin 
lattice models , a doubly degenerate groundstate has been 
detected. Further, in the quantum three-state Potts model 18 , 
a three-fold degenerate groundstate has been detected by us- 
ing a bifurcation of FLS and a probability mass distribution 
function. 

When the system has more than three-fold degenerate 
groundstates, the degenerate groundstates are orthogonal one 
another in the thermodynamic limit. A quantum fidelity de- 
fined by a overlap function between degenerate groundstates 
may not distinguish all the degenerate groundstates properly. 
Then, in this paper, we investigate how to detect generally 
/V-fold degenerate groundstates calculated from a tensor net- 
work algorithm. To do this, we introduce a quantum fidelity 
between degenerate groundstates and an arbitrary reference 
state. The quantum fidelity corresponds to a projection of 
each groundstate onto the chosen reference state. Straightfor- 
wardly, for a broken symmetry phase, the number of different 
projection magnitudes denotes the groundstate degeneracy. At 
a critical point, the different projection magnitudes collapse 
to one projection magnitude. In such a property of the de- 
fined quantum fidelity, the different projection magnitudes of 
the groundstates can be called as multiple bifurcation of the 
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quantum fidelity. A multiple bifurcation point is identified to 
be a critical point. 

As a prototypical example, we explore the groundstate 
wavefunctions in the g-state Potts mode l 19 ' 20 in transverse 
magnetic fields. By employing the infinite matrix prod- 
uct state (iMPS) with the time evolving block decimation 
(iTEBD) method 14 , we calculate the groundstates of the 
model. Due to the broken Z q symmetry, the g-fold degen- 
erate groundstates in the broken symmetry phase are obtained 
by means of the quantum fidelity with q branches. A contin- 
uous (discontinuous) QPT for q < 4 (q > 4) has been mani- 
fested by a continues (discontinuous) fidelity function across 
the critical point. The multiple bifurcation points are shown 
to correspond to the critical points. Also, we discuss a mul- 
tiple bifurcation of local order parameters and its character- 
istic properties for the broken symmetry phase. We find a 
general relation between the order parameters from each of 
degenerate groundstates. The general relation shows clearly 
how the order parameters from each of the degenerate ground- 
states transform under the subgroup of the symmetry group 
Z q . In addition, for q — 4, we calculate the critical expo- 
nents agree well with their exact values. From the finite- 
entanglement scalings of the von Neumann entropy and the 
correlation length, the cental charges are calculated to clas- 
sify the universal classes for each g-state Potts model. 

This paper is organized as follows. In Sec. II, we briefly ex- 
plain the iMPS representation and the iTEBD method in one- 
dimensional quantum lattice systems. In Sec. Ill, the g-state 
Potts model is introduced. Section IV devotes how to detect 
degenerate groundstate by using a quantum fidelity between 
the degenerate groundstates and a reference state. By using 
the quantum fidelity per lattice site, the quantum phase tran- 
sitions are discussed based its multiple bifurcations and mul- 
tiple bifurcation points indicating quantum critical points in 
Sec. V. In Sec. VI, to complete the description of the ordered 
phases, we discuss the magnetizations given from the degen- 
erate groundstates and obtain a general relation between them. 
Section VII presents the critical exponents for q — 4 and the 
central charges for q — 3 and q = 4. Also, we discuss the 
quantum phase transitions from the von Neumann entropies. 
Finally, our summary is given in Sec. VIII. 



II. NUMERICAL METHOD: IMPS APPROACH 

Recently, significant progress has been made in numerical 
studies based on TN representations for the investigation of 
quantum phase transitions— -14 ' 21 , which offers a new perspec- 
tive from quantum entanglement and fidelity, thus providing 
a deeper understanding on characterizing critical phenomena 
in finite and infinite spin lattice systems. Actually, a wave 
function represented in TNs allows to perform the classical 
simulation of quantum many-body systems. Especially, in 
one-dimensional spin systems, a wave function for infinite- 
size lattices can be described by the iMPS representation 14 , 
which has been successfully applied to investigate the prop- 
erties of groundstate wave functions in various infinite spin 
lattice systems 8 ' 11 . 



For an infinite one-dimensional lattice system, a state can 
be written as^ 4 - 



ivp\ _ V V ■ - • jWrK jl'+llrt'+l] i[i+2] . . 



(1) 



where \S ,) denote a basis of the local Hilbert space at the site 
z, the elements of a diagonal matrix A£ are the Schmidt de- 
composition coefficients of the bipartition between the semi- 



infinite chains L(-°o, z) and R(i+ 1, oo), and r e are 
a three-index tensor. The physical indices 5, take the value 
1, • • • , d with the local Hilbert space dimension d at the site i. 
The bond indices ff,- take the value 1, • • • ,x with the truncation 
dimension of the local Hilbert space at the site i. The bond in- 
dices connect the tensors F in the nearest neighbor sites. Such 
a representation in Eq. ([T]) is called the iMPS representation. 
If a system Hamiltonian has a translational invariance, one 
can introduce a translational invariant iMPS representation for 
a state. Practically, for instance, for a two-site translational 
invariance, the state can be reexpressed in terms of only the 
three-index tensors Fa(b) and the two diagonal matrices Aa(b) 
for the even (odd) sites, where {F, A) are in the canonical form, 
i.e., 



\S) [l,r] 



■St-iSiSM---), (2) 



where I and r are the left and right bond indices, respectively. 

Once a random initial state |*F(0)) is prepared in the iMPS 
representation, one may employ the iTEBD algorithm to cal- 
culate a groundstate wavefunction numerically. For instance, 
if a system Hamiltonian is translational invariant and the inter- 
action between spins consists of the nearest-neighbor interac- 
tions, i.e., the Hamiltonian can be expressed by H = > 
where Iv- 1 ' 1 ^ is the nearest-neighbor two-body Hamiltonian 
density, a groundstate wavefunction of the system can be ex- 
pressed in the form in Eq. ©. The imaginary time evolution 
of the prepared initial state 1^(0)), i.e., 



lexpt-HrimO))!!' 



(3) 



leads to a groundstate of the system for a large enough r. By 
using the Suzuki-Trotter decomposition 22 , actually, the imagi- 
nary time evolution operator U = exp[-flr] can be reduced to 
a product of two-site evolution operators U(i, i + 1) that only 
acts on two successive sites z and i+ 1 . For the numerical imag- 
inary time evolution operation, the continuous time evolution 
can be approximately realized by a sequence of the time slice 
evolution gates U(i,i + 1) = exp [-/! [m+1] <5t| for the imagi- 
nary time slice 6t = r/n « 1. A time-slice evolution gate 
operation contracts F^, Fg, one Aa, two Ab, and the evolution 
operator U(i,i + 1) = exp [-/! [,,,+1] c)t|. In order to recover 
the evolved state in the iMPS representation, a singular value 
decomposition (SVD) is performed and the^ largest singular 
values are obtained. From the SVD, the new tensors F A , F B , 
and A A are generated. The latter is used to update the tensors 
Aa as the new one for all other sites. Similar contraction on 
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the new tensors F4, F b , two new Aa, one A B , and the evolution 
operator U(i + 1, i + 2) = exp [-/z [i+1 ' i+2] L and its SVD pro- 
duce the updated F A , F B , and A B for all other sites. After the 
time-slice evolution, then, all the tensors F A , F B , A A , and A B 
are updated. This procedure is repeatedly performed until the 
system energy converges to a groundstate energy that yields 
a groundstate wavefunction in the iMPS representation. The 
normalization of the groundstate wavefunction is guaranteed 
by requiring the norm OFIY) = 1 . Finally, one can determine 
how many groundstates exist for a fixed system parameters 
with different initial states. Once a system undergoes a spon- 
taneous symmetry breaking, the iMPS algorithm can automat- 
ically produce degenerate groundstates with randomly chosen 
initial states for the broken symmetry phase. It has been man- 
ifested by successfully detecting a doubly degenerate ground- 
states from Z2 broken symmetry phases in various spin sys- 
tems such as quantum Ising model, spin- 1/2 XYX with trans- 
verse magnetic field, and a one-dimensional spin model with 
competing two-spin and three-spin interaction ;] 16 ' 17 . However, 
it has not been discussed yet how to detect an A^-fold de- 
generate groundstates. In the following sections, we discuss 
this issue clearly by introducing the g-state Potts model that 
may have a g-fold generate groundstate in broken symmetry 
phases. 
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FIG. 1: (Color online) Groundstate quantum FLS d for the 4-state 
Potts model with an arbitrary reference state as a function of random 
initial state trials. Here, \if>) is an arbitrary reference state, the numer- 
ical groundstates \ifi) are in the broken symmetry phase with A = 0.8, 
and \(p) is in the symmetry phase with A = 1.2, respectively. It is 
clearly shown that, there exist four degenerate groundstates (black 
rectangle) which are labeled by |t/o), \i//2), |<^3), and \ifr4). In the bro- 
ken symmetry phase for A = 1.2, only one groundstate (red circle) 
exists. The dotted lines are guided the eyes. If the number of ran- 
dom initial state trials increases, the probability that the system is in 
each degenerate groundstate approaches 1 /4 in the broken symmetry 
phase for the 4-state Potts model. 



III. QUANTUM (/-STATE POTTS MODEL 

In the lattice statistical mechanics, the quantum g-state 
Potts model is one of the most important models as a gen- 
eralization of the Ising model (q - 2) 20 ' 23 . The g-state Potts 
model has the very intriguing critical behavior that has be- 
come an important testing platform for different numerical 
and analytical methods and approaches in studying critical 
phenomenaSSi 2 ^— . As is well-known in Ref. 20, the quan- 
tum q-state Potts model exhibits a continuous quantum phase 
transition for q < 4 and a first-order (discontinuous) phase 
transition for q > 4 at the critical point. 

We consider the g-state quantum Potts model in a transverse 
magnetic field A on an infinite-size lattice: 



00 q-\ 



M x , p (j) M x , q _ p (j + 1 ) + AM z (j) , (4) 



where A is the transverse magnetic field and M x / z ,p(0 with p e 
[l,q - 1] are the g-state Potts spin matrices at the lattice site 
j. The g-state Potts spin matrices are given as 



A/v, = I 1 7 



and M, = 



q-l 








where I q -\ is the (q - 1) x (q - 1) identity matrix and 
M x ,p - (M Xi {) p . As is known, the model Hamiltonian has a 
Z q symmetry^. If the groundstate of the Hamiltonian does 
not preserve the Z q symmetry, the system undergoes the Z q 
symmetry breaking, i.e., a QPT. Specifically, when the mag- 
netic filed changes across A c = 1, the g-state Potts model un- 
dergoes a QPT between an ordered phase and a disordered 



phase. The QPT is originated from the Z q symmetry breaking 
which results in the emergence of long-range order and the q 
degenerate groundstates in the ordered phase. 



IV. DEGENERATE GROUNDSTATES AND QUANTUM 
FIDELITY 



From a tensor network approach with an infinite lattice 
system, once one obtains a groundstate with the n- 

th random initial state, one can define a quantum fidelity 
F(|i^ (m) ), |0)) = \(if/ (n) \(p)\ between the groundstate and a chosen 
reference state \<p). If F(|i/' ( " ) ), |</>)) has only one constant value 
with random initial states, the system has only one ground- 
state. If F(\i[f (n) ), |0)) has jV projection values with random 
initial states, the system must have degenerate groundstates. 
To distinguish the degenerate groundstates, we employ the 
groundstate FLS in following Ref. |4j as 



lnt/(|(A' 



(n) 



),|0»= lim 

L— >oc 



lnW"),|0)) 



(5) 



where L is the system size. The FLS is well defined in 
the thermodynamic limit even if F becomes trivially zero. 
From the fidelity F, the FLS has several properties as (i) 
d(\if/ {n) ) = \cf>)) = 1 and (ii) its range < d(\^ n) ), \4>)) < 1. 
Within the iMPS approach, the FLS d{\\j/- n) ), \<p)) is given by 
the largest eigenvalue fio of the transfer matrix T up to the cor- 
rections that decay exponentially in the linear system size L. 
Then, for the infinite-size system, d(\i//^" > }, \<p}) = fA " . Actu- 
ally, in order to study quantum critical phenomena in quantum 
lattice systems, Zhou and Barjaktarevic defined the FLS from 
the fidelity between groundstates^^. The Zhou-Barjaktarevic 
FLS has been manifested as a model-independent and uni- 
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FIG. 2: (Color online) Groundstate quantum FLS d for the 4-state 
Potts model as a function of random initial state trials. The reference 
states are chosen in (a) |Ti>, (b) |T 2 >, (c) |¥ 3 >, (d) |^ 4 >, and (e) \<p). 
It is clearly shown that, if the reference state is chosen as one of 
the degenerate groundstates in the broken symmetry phase, the other 
three groundstates has the same value of the FLS in (a)-(d). If the 
reference state is chosen as the groundstate in the symmetry phase, 
the FLSs for the four degenerate groundstates has the same value in 
the broken symmetry phase in (e). 



versal indicator successfully detecting quantum phase tran- 
sition points including KT transition and topological phase 
transition 11 . The FLS introduced in Eq. (O for this study is 
a simple extension of the Zhou-Barjaktarevic FLS. By means 
of a bifurcation of the Zhou-Barjaktarevic FLS and a probabil- 
ity mass distribution function, one can distinguish degenerate 
groundstates and detect a quantum phase transition. 

In order to show clearly how to detect degenerate ground- 
states based on Eq. (0, as an example, we consider the 4-state 
Potts model. In Fig. [1] we plot the FLS d for random initial 
states. Here, A = 0.8 and 1.2 are chosen. Due to the broken Z4 
symmetry for q — 4, the system has 4 degenerate groundstats 
for the broken symmetry phase A < Ac = 1. It is clearly shown 
that, for A = 0.8 (A < A c ), there exist four different values of 
the FLS, while, for A = 1.2 (A > A c ), there exists only one 
value of the FLS. From each value of the FLS for A = 0.8, 
we label the four degenerate groundstates as \i/fi), |ife)> I^X 
and It/14). For for A = 1.2, the groundstate is denoted by \ip). 
Actually, we have chosen more than 200 random initial states. 
The probability P(n), that the system is in each groundstate 



for A < A c , is found to be P(n) 1 /4 in the broken symmetry 
phase. Then, for the ^-state Potts model, with a large number 
of random initial state trials, one may detect the q degener- 
ate groundstates with the probability P{n — > 00) = l/q finding 
each degenerate groundstate in the broken Z q symmetry phase. 

One may choose the reference state as one of the degen- 
erate groundstates. In Fig. [2j we choose the reference state 
as one of the groundstates for the broken symmetry phase in 
(a) l^i ), (b) \W 2 ), (c) |T 3 ), and (d) PF 4 >, and as the ground- 
state \<p) for the symmetry phase in (e). In Figs. HJa)-(d), 
it is clearly shown that, for A = 0.8 (A < A c ), there exist 
two different values of the FLS, i.e., d(\i/r m ), \i/f m )) = 1 and 
d(\ifr m ), l^m^m')) = constant + 1 (m,m' e [1,^]). Further- 
more, d(\ifr m ), \(p)) = constant. In Fig. [2] (e), the reference 
state chosen as the groundstate in the symmetry phase shows 
d(\ip), \(p)) = 1 and d(\ifr m ), \(p)) - constant + 1. In this case, as 
discussed in Ref. [jjl, one may use a probability mass dis- 
tribution function, as an alternative way, to distinguish the 
degenerate groundstates. Here, it is shown that the degener- 
ate states can be distinguished by a reference state chosen as 
arbitrary states except for the degenerate groundstates in the 
broken symmetry phase and the groundstate in the symmetry 
phase. 



V. QUANTUM FIDELITY PER LATTICE SITE FOR 
PHASE TRANSITIONS 

In the view of a groundstate degeneracy, the degenerate 
groundstates in the broken symmetry phase exist until the sys- 
tem reaches its critical point. Once one can detect degener- 
ate groundstates, one may also detect a quantum phase transi- 
tion by the quantum fidelity in Eq. (f5]). Detected degenerate 
groundstate wavefunctions may allow us also to investigate di- 
rectly a property of quantum phases even for unexplored bro- 
ken symmetry phases. In this section, then, we will discuss 
the quantum phase transitions for the g-state Potts model. In 
the next section, we will discuss the relation between system 
symmetry and order parameter directly from groundstates. 

From the detected degenerate groundstates, in Fig. [3] we 
plot the FLSs as a function of the transverse magnetic field A 
for q = 3, 4, and 5. Here, the truncation dimension is chosen 
a&x - 32. Figure [3] shows clearly that, as the transverse mag- 
netic field decreases, the FLSs in the symmetry phase branch 
off the q FLSs. The branch points are estimated numerically 
as A c = 1.0004 for the 3- and 4-state Potts models. Also, for 
the 5-state Potts model, the multiple bifurcation point exists 
exactly at A c = 1 . These results show that the branch points 
agree well the exact critical point A c — 1 . Then, the branch 
points correspond to the critical point. Such a branching be- 
havior of the FLS can be called as multiple bifurcation and a 
branch point as multiple bifurcation point. Consequently, it is 
shown that the FLS from the quantum fidelity between degen- 
erate groundstates and a reference state can detect a quantum 
phase transition. 

As is known, for the g-state Potts model, the continues (dis- 
continuous) phase transitions occurs for q < 4 (q > 4). Figs. 
[3](a) for q = 3 and (b) q = 4 show that the FLSs are a con- 
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FIG. 3: (Color online) Groundstate quantum FLS ii for (a) 3-, (b) 
4-, and (c) 5-state Potts models as a function of the transverse mag- 
netic field A with the truncation dimension x = 32. In the broken 
symmetry phase, the q branches of the FLS correspond to the num- 
bers of the q degenerate groundstates. As the magnetic field crosses 
the critical point A c , the FLSs show the multiple bifurcations with 
three, four, and five branches, respectively, in the broken symmetry 
phase. Note that, for the 3- and 4-state Potts models, the FLSs show 
their continues behaviors, while, for the 5-state Potts model, the FLS 
shows a discontinues behavior. These continues and discontinues be- 
haviors indicate the continues phase transitions for q = 3 and 4, and 
the discontinues phase transition for q = 5. The numerical multiple 
bifurcations are given as A c = 1.0004 for q = 3 and 4 and A c = 1 for 
q = 5. 



tinues function for the quantum phase transition. While Fig. 
[3](c) q = 5 shows that the FLS is a discontinues function for 
the quantum phase transition. Then, the continues (discon- 
tinues) behaviors at the critical points imply that a continues 
(discontinues) phase transition occurs. As a result, the FLS 
in Eq. Q can clarify continues (discontinues) quantum phase 
transitions. 
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FIG. 4: (Color online) Magnetization (a) (M v ,) and (b) {M Xt2 ) as a 
function of the transverse magnetic field A for the 3-state Potts model. 
For the broken symmetry phase A < A c , each magnetization is given 
for each of the three degenerate groundstates denoted by \if/ m ) (m £ 
[l,q]. The numerical critical point locates at A c = 1.0004. 



VI. SYMMETRY AND ORDER PARAMETER 

The existence of a degenerate groundstate means that each 
of degenerate groundstates possess its own order described by 
each corresponding order parameter. The each order parame- 
ter, which is nonzero value only in an ordered phase, should 
distinguish an ordered phase from disordered phases—. To 
complete the description of an ordered phase, it is required 
to specify how each order parameter from each of degener- 
ate groundstates transforms under a symmetry group G that 
is possessed by the Hamiltonian because each order parame- 
ter is invariant under only a subgroup of the symmetry group 
G although the Hamiltonian remains invariant under the full 
symmetry group G 2 . Then, in this section, we discuss local 
magnetizations obtained from each of the degenerate ground- 
states. 

Let us first discuss the local magnetizations for the 3-state 
Potts model, i.e., q = 3. In Fig. |4] we plot the magneti- 
zation (a) (M x \) and (b) {M x t) as a function of the traverse 
magnetic field A. The magnetizations disappear to zero at the 
numerical critical point A c = 1 .0004. Also, all the magneti- 
zations show that the phase transition is a continuous phase 
transition. For the broken symmetry phase A < A c , each 
magnetization is calculated from each of the three degenerate 
groundstates, which is denoted by \tf/ m ). Note that all the ab- 
solute values of the magnetizations (ifr m \M x<p \iJf m ) = (M xp ) m 
are the same values at a given magnetic field. Furthermore, 
for a given magnetic field, the magnetizations in the com- 
plex magnetization plane seem to have a relation between 
them under a rotation, which is characterized by the value 
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FIG. 5: (Color online) Magnetizations (a) (M Xj i), (b) (M v _2>, and (c) (M x ^) as a function of the transverse magnetic field A for the 4-state 
Potts model. For the broken symmetry phase, the magnetizations are given from each of the four degenerate groundstates. The critical point is 
estimated numerically as A c = 1.0004. 



a>3 = exp[2m/3]. Then, in Figs. Ufa) and (b), it is observed 
that, for a given magnetic field A < A c , there are the rela- 
tions between the magnetizations as {M x ^)\ = a>j 1 (M^i )2 = 

ojf(M Xt ih and <M. t>2 )i = ^ 2 (M xa h = ^(M^h- Also, 
for each groundstate wavefunction, the magnetizations have 
the relations as (M.u)i = (M x , 2 )u {M xA ) 2 = oi^{M x ,2)i, 
and (M A j) 3 = o/T 2 \M Ai2 ) 3 . These results show that, in the 
complex magnetization plane, the rotations between the mag- 



Next, we consider the magnetizations for the 4-state Potts 
model, i.e., q — 4. In Fig. [5] we plot the magnetization (a) 
(M x ,i), (b) (M A>2 ), and (c) (M Aj3 } as a function of the traverse 
magnetic field A. The numerical critical point is A c = 1.0004. 
All the magnetizations show that the phase transition is a con- 
tinuous phase transition. Similar to the case q — 3, all the 
absolute values of the magnetizations have the same values at 
a given magnetic field and the magnetizations in the complex 
magnetization plane have a relation between them under a ro- 
tation, which is characterized by the value a>4 = exp[2m/4]. 
In Fig. [5] we observe that, for a given magnetic field A < 
A c , there are the relations between the magnetizations as 
(M xA h = ^{M^h = w 4 2 (M.u> 3 = « 4 3 <M. U ) 4 from (a), 
(M xa h = ajf{M xa ) 2 = u 4 4 {M xa h = w 4 6 <M A , 2 >4 from (b), 
and (M x , 3 h = w 4 3 <M A>3 > 2 = w 4 6 (M x , 3 ) 3 = w 4 9 (M A -, 3 ) 4 from 
(c). Also, for each groundstate wavefunction, the magneti- 
zations have the relations as {M x> \)\ = (M t>2 )i - {M x ^)\, 
(M xA ) 2 = oj- x {M x . 2 ) 2 = ojf{M x , 3 ) 2 , (M X>1 ) 3 = cjf(M xa h = 
w 4 4 <M Xj3 ) 3 , and (M xA } 4 = w 4 3 <M Aj2 ) 4 = w 4 6 (M Aj3 ) 4 . These 
results show that, in the complex magnetization plane, the 
rotations between the magnetizations for a given magnetic 



netizations for a given magnetic field are determined by the 
characteristic rotation angles 0-0, 2n/3, and 47r/3, i.e., 
(M x , p ) m = g 3 {M XtP ,) m , with ^ 3 6 {/,w 3 , w 2 }. Here, we have 
chosen the \tff\) that gives a real value of the magnetizations, 
i.e., (M x i)i and (M x 2 )\ are real. Also, the degenerate ground- 
states give the same values for the z-component magnetiza- 
tions, i.e., (M z )i = (M z ) 2 = (M_> 3 . 



I 

field are determined by the characteristic rotation angles = 
0, 27r/4, 47r/4, and 67r/4, i.e., (M XtP ) m = g4(M x ^) m , with 
g4 e {I, a>4, a» 4 , w 4 }. Also, the degenerate groundstates give 
the same values for the z-component magnetizations, i.e., 
(M z h = {M : ) 2 = (M z ) 3 = (M z ) 4 . 

From the discussions for the 3-state and 4-state Potts mod- 
els, one may refer a general relation between the magnetiza- 
tions for any g-state Potts model. Actually, for any q, we find 
the relations 

af g a ~ m) (M XtP ) m = uf x - m '\M x , p ,) m ,, (6a) 
(M z ) m = (M : ) m ,, (6b) 

where a) q = exp[27n'/g]. When one calculate the magnetiza- 
tions of the operator M xp with different wavefunctions, the 
relations in Eq. (O reduce to {M xp ) m = oj^ m m \M XtP ) m i 
that satisfies the relations of the magnetizations in Figs. |4]and 
Also, if one choose a wavefunction, the magnetizations of 
the operators M x j, ■ ■ ■ , M Ai9 _ 2 , and M x>q -\ have the relations 
{M xp ) m = idq p)i ~ l m) {M x<p >) m reduced from the relations in 
Eq. (0. Furthermore, these results show that, in the complex 
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FIG. 6: (Color online) Magnetizations (a) (M v j), (b) (M x o), (c) (M x ^), and (d) (M v _4> as a function of the transverse magnetic field A for 
the 5-state Potts model. For the broken symmetry phase, the magnetizations are given from each of the five degenerate groundstates. The 
numerical critical point is the exact value A c = 1 . 



magnetization plane, the rotations between the magnetizations 
for a given magnetic field are determined by the characteristic 
rotation angles = 0, 2n/q, 4-n/q, 6n/q, • ■ • , and 2(q - \)nlq. 
Thus, Eq. (|6]l can be rewritten as 



(M XiP ) m = g q (M x , p ,)„ 

,2 



(7a) 

{/.Wj.i^-.-.wf 1 }. (7b) 

Then, Eq. (|7]l shows clearly that the g-state Potts model 
has the discrete symmetry group Z q consisting of q elements. 
For q — 2, the g-state Potts model becomes the quantum 
Ising model that has a doubly degenerate groundstate due 
to a Z2 symmetry. From Eq. 0, one can easily confirm 
(M,)i = -{M x )2 because of g2 e {/, coi) with a>2 = exp[m']. 
Also, the g-state Potts Hamiltonian in Eq. (HJi is invariant with 
respect to the q way transformations, i.e., U, n H q uj n = H q for 
m e [l,q], which implies that the system has the g-fold degen- 
erate groundstates for the broken symmetry phase according 
to the spontaneous symmetry breaking mechanism. The trans- 
formations are given as 



M 



M 7 



4) m_1 - 
M, 



(8) 



where a> q = exp[2ni/q]. Although the g-state Potts Hamil- 
tonian remains invariant under the full q transformations, for 
the Z q broken symmetry phase, the q order states described 
by the q degenerate groundstates are invariant under only the 
subgroup of the Z q symmetry group. Obviously, Eqs. © and 
(0 show the relations between the order parameters of the q 
equivalent ordered states under the q transformations. As a 
consequence, it is shown that, from the degenerate ground- 
states, one can determine the order parameters as the magne- 
tizations and their specification of how the order parameters 
transform under the symmetry group Z q . 



To make clearer the general relation of the magnetizations 
in Eqs. (O and ©, let us consider the 5-state Potts model. 
In Fig. [6] we plot the magnetization (a) {M Xt \), (b) (M x o), 
(c) {M x ,t), and (d) (M v .4) as a function of the traverse mag- 
netic field A. The numerical critical point is obtained as the 
exact value Ac = 1. Also, all the magnetizations shows that 
the phase transition is a discontinues phase transition. Simi- 
lar to the cases q — 3 and q = 4, all the magnetizations have 
the same values at a given magnetic field, and the magne- 
tizations in the complex magnetization plane have a relation 
between them under a rotation, which is characterized by the 
value u>5 = exp[2m'/5]. From Figs. |6](a)-(d), the relations 
between the magnetizations agree with those in Eqs. (0 and 
0. 



VII. CRITICAL EXPONENTS, CENTRAL CHARGE, AND 
UNIVERSALITY 



As is known, for q < 4, the quantum phase transitions 
are a continuous phase transition. If q > 4, discontinuous 



TABLE I: Critical exponents and central charge for 4-state Potts 
model 



Exponents 


a 


P 


7 


6 


V 


1 


c 


Exact" 
iMPS 


2/3 


1/12 
0.0843 


7/6 
1.0718 


15 


2/3 
0.6300 


1/4 
0.2510 


1 

0.9803 



"The exact vaules of the critical exponents and central charge are taken from 
Refs.HJandlll 
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FIG. 7: (Color online) (a) Magnetization (M x ) and (b) susceptibility 
d,i(M x ) as a function of \A — A c \ near the critical point A c = 1.0004 for 
4-state Potts model. 



(first order) phase transitions occur in the g-state Potts model. 
In this sense, then, the critical q is q c = 4. In the section, 
we will study the critical exponents of q c — 4 based on the 
iMPS groundstate wavefunctions. Actually, all the degenerate 
groundstates for the broken symmetry phase give a same ex- 
ponents as it should be. In addition, in order to calculate a cen- 
tral charge, we discuss the von Neumann entropy for q — 2>,A, 
and 5. In the Table U for the 4-state Potts model, the critical 
exponents and the critical charge from the our iMPS results 
are compared with their exact values. 

Critical exponents - In our iMPS approach, we obtain the 
4-fold degenerate groundstates at zero temperature. The criti- 
cal exponents a for specific heat and 5 for the field dependence 
of the magnetization at the critical temperature cannot be cal- 
culated. In the following, we discuss the four exponents, i.e., 
P, y, rj, and v. Thus, we start with the magnetization near the 
critical point A c = 1.0004. In Fig. [7] (a) the magnetization 
(M x ) and (b) the susceptibility d^(M x ) are plotted as a func- 
tion of \A - A c \ in the log-log plot. It is shown clearly that both 
the magnetization and the susceptibility can be described by 
power laws, i.e., (M x ) cc \A-A c f and d A (M x ) oc \A-A c \- y with 
their characteristic exponents p and y, respectively. Thus, the 
fitting function is chosen to be \n(M x ) = pin\A - A c \ + /3o 
for the magnetization and ]nd^(M x ) = -yln \A - A c \ + yo for 
the susceptibility with the fitting constants /3o and yo, respec- 
tively. From the numerical fittings, we obtain (a) p = 0.0843 
and Po = 0.0394 for the magnetization and (b) y = 1.0718 
and yo = -3.0308 for the susceptibility. The fitted critical 
exponents p = 0.0843 and y = 1.0718 are quite close to the 
exact values 20 p = 1/12 (= 0.0833) and y = 7/6 (= 1.1667), 
respectively. 

In Fig. [8] we plot (a) the correlation length £ as a function 
of \A— A c \ and (b) the correlation (M x (i)M x (j)) as a function of 
the lattice site distance at the critical point A c = 1 .0004. It 
is shown clearly that both the correlation length and the corre- 
lation at the critical point A c can be described by power laws, 
i.e., £ oc \A - A c \- y and {M x (i)M x (j)) oc\A- A c \-' ! with their 
characteristic exponents v and 77, respectively. We choose the 
fitting functions as ln£ = -vln \A — Ac\ + vo for the correlation 
length and ln(M x (i)M x (j)) = —77 In |i - j\ + 770 for the cor- 



FIG. 8: (Color online) (a) Correlation length £ as a function of \A - 
A c \ for the 4-state Potts model, (b) Correlation (M x (i)M x (j)) as a 
function of the lattice site distance \i - j\ at the critical point A c = 
1.0004. 



relation at the critical point with the fitting constants vo and 
770, respectively. From the numerical fittings, we obtain (a) 
v = 0.6300 and vo = -1.2153 for the correlation length and 
(b) 77 = 0.2510 and 770 = -3.5644 for the correlation. The 
fitted critical exponents v = 0.6300 and 77 = 0.2510 are quite 
close to the exact values 20 v — 2/3 and 77 = 1/4, respectively. 

von Neumann entropy and central charge.— In our iMPS 
approach, the von Neumann entropy S can be directly eval- 
uated by the elements of the diagonal matrix A„. that are 
the Schmidt decomposition coefficients of the bipartition be- 
tween the semi-infinite chains L(-oo, 7) and R(i + 1, 00). 
This implies that Eq. (fl]i can be rewritten by ^ = 
H«=i ^a#»>l<Aff}, where and |t/>*> are the Schmidt bases 
for the semi-infinite chains L(~co, /) and R(i + 1, 00), re- 
spectively. For the bipartition, then, the von Neumann entropy 
S can be defined as ?1 S — -Tt[ql logg/J = -Tr[p« logp«], 
where qi = Tr« q and qr = Tr^ q are the reduced density ma- 
trices of the subsystems L and R, respectively, with the den- 
sity matrix g = | < P)( V F|. For the semi-infinite chains L and 
R in the iMPS representation, the von Neumann entropy S 
calculated by S = - _\ A%logA%. Actually, the von Neu- 
mann entropy as one of the quantum entanglement measures 
have been proposed as a general indicator to determine and 
characterize quantum phase transitions 32,33 . Also, the loga- 
rithmic scaling of the von Neumann entropy was conformed 
to exhibit conformal invariance and the scaling is governed 
by a universal facto r 5,34 ' 35 , i.e., a central charge c of the as- 
sociated conformal field theory. In the iMPS approach, for a 
continuous phase transition, the diverging entanglement at a 
quantum critical point gives simple scaling relations 34 for (i) 
the von Neumann entropy S and (ii) a correlation length f v 
with respect to the truncation dimension^ as S ~ log^ and 
f v ~ Ax", where c is a central charge, k is a so-called finite- 
entanglement scaling exponent, and A is a constant. By using 
the relations, a central charge can be obtained numerically at 
a critical point. 

In Fig. |9](a), we plot the von Neumann entropies as a func- 
tion of the transverse magnetic field A for the 3-, 4-, and 5-state 
Potts models with the truncation dimension ^ = 32. In the en- 
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FIG. 9: (Color online) (a) Groundstate von Neumann entropies for 
the 3-, 4-, and 5-state Potts models as a function of the transverse 
magnetic field A. (b) Von Neumann entropies 5 and (c) Correlation 
lengthes f v for q = 3 and q = 4 with respect to the truncation dimen- 
sion ^ at the exact critical point A c = 1 . 



tropies, there are singular points at A c = 1 .0004 for g = 3 and 
4, and A c — 1 for g = 5, which are consistent with the criti- 
cal points from the multiple bifurcation points in Fig. [3] and 
from the magnetizations as the order parameters in Figs. |H 
[5] and [6] It is shown clearly that the von Neumann entropies 
capture the phase transitions. The (dis-)continuity of the von 
Neumann entropies for q — 3 and q — 4 (q — 5) indicates a 
(dis-)continuous phase transition between the broken symme- 
try phases and the symmetry phases. However, the von Neu- 
mann entropies for the different degenerate groundstates give 
the same values, which implies that the von Neumann entropy 
cannot distinguish the different degenerate groundstates in the 
broken symmetry phases. 

For q = 3 and q = 4, in Figs. [9] (b) and (c), we plot 
(b) the von Neumann entropy and (c) the correlation length 
f v as a function of the truncation dimension x at the critical 
points A c = 1. Here, the truncation dimensions are taken 
as x = 8 ' 16, 25, 32, 50, and 64. It is shown that both the 



von Neumann entropy S and the correlation length £ v diverge 
as the truncation dimension x increases. In order to obtain 
the central charges, we use the numerical fitting functions, 



i.e., S q (x) 



+ b q \ogx and &,q(x) = A qX Kq - Numeri- 



cally, the constants of the von Neumann entropies are fitted 
as a 3 = -0.1215 and£ 3 = 0.2000 for q = 3 anda 4 = -0.1830 
and Z?4 = 0.2016 for q = 4. The power-law fittings on the cor- 
relation lengthes g v$q give the fitting constants as A3 = 0.2060 
and k 3 = 1.4775 for q = 3. and A 4 = 0.2243 and k 4 = 1.2340 
for q = 4. As a result, from the K q and b q , the central charges 
are obtained as C3 = 0.8121 for q — 3 and c\ = 0.9803 for 
q — 4, which are quite close to the exact values 30 c = 4/5 for 
q — 3 and c = 1 for q = 4, respectively. 

VIII. SUMMARY 

We have investigated how to obtain degenerate ground- 
states by using a quantum fidelity. To do this, we have in- 
troduced the quantum fidelity between the degenerate ground- 
states and an arbitrary reference state. The distinguished de- 
generate groundstate wavefunctions allows naturally to detect 
a quantum phase transition that is indicated by a multiple bi- 
furcation point of the quantum fidelity per lattice site. Fur- 
thermore, as the complete description of an order phase, it is 
possible to classify how the order parameters calculated from 
the degenerate groundstates transform under the subgroup of 
a symmetry group of Hamiltonian. As an example, the g-state 
Potts model has been investigated by employing the iMPS 
with the iTEBD. We have obtained the g-fold degenerated 
groundstates explicitly. A general relation between the mag- 
netizations calculated from the degenerate groundstates is ob- 
tained to show the spontaneous symmetry breaking of the Z q 
symmetry group in the g-state Potts model. In addition, the 
critical exponents and the central charges are directly calcu- 
lated from the degenerate groundstates, which is shown that 
the iMPS results are quite close to the exact values. 
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